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Abstract 

This paper is devoted to the investigation of modified Chaplygin gas model 
in the context of solvable k-essence cosmologies. For this purpose, we con- 
struct equations of state parameter of this model for some particular values of 
the parameter n. The graphical behavior of these equations are also discussed 
by using power law form of potential. The relationship between k-essence 
and modified Chaplygin gas model shows viable results in the dark energy 
scenario. We conclude that the universe behaves as a cosmological constant, 
quintessence phase or phantom phase depending upon n. 



1 Introduction 

The accelerating cosmic expansion of the universe, motivated by the mysterious 
dark energy (DE) idea, was noticed in 1998 through supernovae experiments [1]. A 
number of various cosmological observations claim that this energy fills the space 
but yet there is no clue about its identity. The expansion of the universe started 
to accelerate when the energy density of matter was overcome by the DE. The 
evolution of DE density is governed by an equation of state parameter, uj = p/ p. This 
parameter has raised recently many interesting features for accelerating expansion 
of the universe. Some of them are as follows: 

• Qunitessence region obeys — l<Ci;<— 1/3[2], 
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Cosmological constant is considered to be the simplest case for DE for which 

u = —1, 



• For phantom phase, the EoS parameter follows the domain u < —1 |3]. 

• Qunitum model inherits both of the properties of quintessence and phantom 
phases by crossing the phantom divide line u = —1 

Modified Chaplygin Gas (MCG) model [5j is one of the proposed models to 
discuss the DE dominated phase for accelerating expansion of the universe. It is 
a single fluid model to unify the DE and dark matter. This model behaves as 
dark matter when its energy density evolves as p oc in the early times of the 
universe (where the scale factor a is a function of time), whereas the constant energy 
density behaves as DE era in late time acceleration. Thus the MCG model provides 
unification between dark matter and DE |6j. The scalar field model is another 
category to discuss the DE scenario. Recently, k-essence [7j gains an important role 
to explore the djTiamics of the universe. It inherits a kinetic term which deals with 
late time inflation and early time acceleration named as k- inflation [8J. 

Another DE model, g-essence, is a general essence model reduced in two models: 
k-essence (kinetic quintessence) and f-essence (fermionic k-essence) [5]. Razina et 
al. [10] explored the relationship between fermions and scalar fields by taking two g- 
essence models. They concluded that g-essence describes the accelerating expansion 
of the universe. Tsyba et al. [11] reconstructed the interaction among Chaplygin 
gas. Generalized Chaplygin gas as well as fermions and found EoS parameter for 
particular f-essence models. Myrzakulov [12] discussed the fermionic k-essence to 
explore accelerating expansion of the universe. Jamil et al. [13] investigated the 
behavior of solvable f-essence cosmology in MCG scenario. 

In this paper, we study the purely kinetic k-essence cosmology by using the 
phenomenon of MCG model. We obtain a non-linear differential equation which 
depends upon the parameter n. The EoS parameter is evaluated by using some 
particular values of n and is displayed graphically to investigate the dynamics of 
the universe. The power law form of potential function is used to obtain the kinetic 
term of purely kinetic k-essence. 

The paper is organized as follows. In the next section, we present some basics of 
k-essence model and the MCG model. In section 3, the EoS parameter is constructed 
for different cases and their graphical representation is discussed. Section 4 describes 
the results for some other forms of solvable k-essence cosmology. In the last section, 
we summarize and conclude the results. 

2 K-essence and MCG Formalism 

Here, we consider the class of models described by the following action 




(2.1) 



where is the scalar field, R is the Ricci scalar. 
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is the canonical kinetic energy and gf^i, is the metric tensor, is the covariant 
derivative associated with metric g^u- We are interested in flat, homogeneous and 
isotropic FRW universe described by the line element 

ds^ = -dt^ + a^{dx^ + dy"^ + dz^), (2.2) 

where a{t) is the scale factor. The k-essence equations corresponding to the action 
(12. ip read as 

3H^-p = 0, (2.3) 
2H + 3H^+p = 0, (2.4) 
Kx4>+iKx + 3HKx)^-K^ = 0, (2.5) 
p + 3H{p + p) = 0. (2.6) 

Here dot denotes derivative with respect to cosmic time, the Hubble rate is H = d/a, 
the kinetic term, the energy density and the pressure, respectively, are 

X = O.50^ p = 2XKx-K, p = K. (2.7) 



It was shown in [13] that for this model the speed of propagation of scalar pertur- 
bations (speed of sound), c^, is given by 

px PX + 2XPX ^ ' 

In the case of purely kinetic k-essence, K = K{X), Eg. (12. 5 p takes the form 

Kx4>+{Kx + 3HKx)^ = 0, 

yielding the general solution 

XKl = ka-\ (2.9) 

with A; > [ISIIIS]. 

The EoS of the MCG DE model is [IT] 

p = Ap-^, (2.10) 

where A and B are positive constants and < a < 1. Using Eqs.f l2.6p and (I2.10p . 
the MCG energy density and pressure evolve as 

p = [5(l + v4)-^ + Ca-=^(i+")(i+^)]^, (2.11) 
p = [ACa-=^(^+")(^+^) - 5(1 + A)-'] [B{1 + A)'' + Ca-=^(^+")(^+^)] "^2.12) 
In this case, the EoS parameter is 

_ ACa-3(i+")(i+-^) - B{1 + Ay^ 

5(1 + A)-i + Ca-3(i+a)(i+A) ■ ^2.13) 

Here we assume the potential in the form of power law [18] 

= U^. (2.14) 
Using the above equation, the kinetic term takes the form 

X = 0.502/3 V(^-i). (2.15) 
For the sake of simplicity, we assume 00 = 1. 
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3 Solvable k-essence cosmology: K = F[X) 

In this section, we would like to consider the case [12] when 

K = F{X) (3.1) 

that is the purely kinetic k-essence case. We solve the following system 

p = 2XFx~F, (3.2) 

F = Ap-^. (3.3) 

Solving Eq.(IS2]) for F, it follows that 

F = EVx + y^l^^dX, (3.4) 

where E is an integration constant. Equations (13. 3 p and (13. 4p give 



eVx + ^ [^dX = Ap-^ 
2 J p^ 



which has the solution 



(1 + A)pi+° - (W^X)'^p"(^+") -5 = 0, (3.5) 
where is a constant and 

n-^ii±^ (36) 

From Eq. (l3.6p . it follows that A and a are related by 

in — l)a nA 
^ = ~i 7^ or a= — --. (3.7) 

{n-l)a + n [n - + A) ^ ' 

The search for the analytical solutions of Eq. (l3.5p is a hard job. We would like 
to find particular solutions for some values of n to get insights of the model. It is 
mentioned here that the values of constants [131 US] in the models are taken in such 
a way that we obtain an expanding universe with acceleration. 

3.0.1 Example 1: = 

First, we consider the particular value n = 0. It follows from Eq. (l3.6p that this value 
yields two cases, either a; = 0ory4 = — 1. 

Case 1 

Let us consider the case a = 0. Consequently, Eqs. (13.20 and (13. 3p take the form 

p = 2XFx -F, F = Ap~ B. (3.8) 
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Figure 1: Plot of tu versus t with A = —0.1, W = 0.2, B = 0.5. In the left graph, 
/3 = —2 (black), (3 = —4 (red), f3 = —6 (green), /3 = — 10 (blue) and in right graph 
(3 = 2 (black), /3 = 4 (red), (3 = 6 (green), (3 = 10 (blue). 



Solving the above equations, the equation for p reads as 

(1 + A)p - (WX)^ - £ = 0, 



which yields 



p= {1 + A)-'^[B + (WX)"^]. 
The expression for pressure becomes 



P 



1 + A)-'[-B + A{WX)—]. 



The corresponding EoS parameter is given by 

B{l + A) 



u = A- 



B + iWXy-i^ 



Using Eq. fl2.15p . the above equation can be written as 



A~ 



B{l + A) 



5+ (iy/32t2(/3-l)/2)^' 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



The behavior of EoS parameter uj in terms of time is shown in Figure 1 which 
indicates the cosmological evolution. For /3 > 0, initially uj shows the DE phase of the 
universe dominated by quintessence. As the value of (3 increases, the EoS parameter 
becomes constant aX u = —1 indicating the universe behavior like cosmological 
constant. On the other hand, for /3 < 0, the EoS parameter possesses negative 
values and becomes constant at —0.1. It presents a universe which just evolves from 
matter dominated universe [u = 0) to DE energy dominating universe. 

Case 2 



Now we consider the case when A = —1. The corresponding equations (13.21) and 
(13. 3p read as 

p = 2XFx -F, F = -p- Bp-"", (3.13) 
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Figure 2: Plot of u versus t with W = 0.2 and /3 < 0. 



which lead to 

aB In p - {l + a)-^p^+'^ + \n{WX)^ = 0. (3.14) 

Consider some particular solutions of this equation. 

i) If a = 0, then we have 

p = \n{WX)^ (3.15) 

and for the pressure 

p = \n{WX)-^ - B. (3.16) 
The corresponding EoS parameter in terms of potential is given by 

2 

The EoS parameter u shows the radiation and matter dominated universe for t < 1.3 
as shown in Figure 2. After this interval, the universe converges to DE dominated 
phase lying in the quintessence region as— l<a;<— 1/3. 

ii) For a = —1, the energy density and pressure become 

p={WX)^(^, p= -{1 + B){WX)^oTB). (3.18) 

The EoS parameter takes the following form 

uj = -l-B. (3.19) 

For this case, the behavior of the universe depends upon the value of constant B. 
If we take B = 0, the universe converges to DE phase dominated by cosmological 
constant. For B > 0, the universe meets the phantom region while B < shows a 
quintessence DE universe. 

3.0.2 Example 2: n = 1 

If n = 1, Eq. fl3.6p yields A = 0. This case corresponds to the generalized Chaplygin 
gas model [16\ for which the expressions of p and p are as follows 

p = B^[l-{WX)^]^, (3.20) 
p= -B^[l - (WX)^]^. (3.21) 



6 



(0 



-0 . 985 
-0 . 9875 
-0 . 99 
-0 . 9925 

-0 . 995 
-0 . 9975 




Figure 3: Plot of to versus t with W = 0.2, (3 = -2 and a = 0.2 (black), a = 0.4 
(red), a = 0.6 (green). 



The corresponding EoS parameter is given by 

uj = -i + (WX)^. 
Inserting the value of X, Eq.f l3.22p becomes 

OJ = —1 + 



(3.22) 



(3.23) 



Figure 3 shows the graphical behavior of the EoS parameter with respect to evolving 
time. The universe acquires quintessence region initially. As the time elapses, it 
becomes constant at u = —1 showing a universe dominated by the cosmological 
constant. 

3.0.3 Example 3: n = 2 

Now we consider the case when n = 2. Here A and a are related by the formula 



A 



a 



2A 



or a 



l + A 



a + 2 

The equation for p (13.51) takes the form 

(1 + A)pl+" - (VrX)^p2(l+a) _ 5 = 0, 

which has the solution 



p = {wxy 



^Tl\/ii + Ay-ABiwxy 



or 



p = {WX)-T. 
The pressure is given by 

a{WX)- 



2 + a 



It yi- B{2 + ay{wxy 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



P 



a + 2 [2 + a 



1 



iTyl- B{2 + ay{wxy 
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Figure 4: Plot of uj versus t with W = 0.2, B = -0.1, a = 1.5 and /3 < 0. 



BWX 



1 



2 + a 

The corresponding EoS parameter is 



u = ^ - B{WX)^ 



It yi- B{2 + ay{wxy 



(3.28) 



1 + Q 



1 ^ -y/l - B{2 + ay{WX)^ 



-1 



a + 2 [l + a 

Inserting the value of potential and after simplifying the positive root, it yields 

2S(l + a)(iy/32t2(/3-i))^ 



(3.29) 



a 



a + 2 2- 



B{a + 2y{WPH^(^~ 



1)^ 



(3.30) 



For this case, the EoS parameter results a universe having dynamics of quintessence 
of DE scenario shown in Figure 4. It becomes constant at = —0.4285 with the 
passage of time. The negative root of Eq. fl3.29p does not give any fruitful result. 



3.0.4 Example 4: n = 0.5 

If n = 0.5 then A and a satisfy the relation 

a 



A 



A 



or a 



a-1 " " l+A 
The equation (13. 5 p becomes 

(1 + A)pi+" - (l^X)^p°-^(^+") -5 = 0. 

The solution of this equation reads as 



P=[2{l + A)\ 



l + a 



{WX)— ± y {WX)^ +AB{1 + A) 



or 



l-a 



2 

l + a 



{WX) 



l + g 
Aa T. 



iWX)^ + 



AB 
1 — a 



2 

l + a 



2 

l + a 



(3.31) 



(3.32) 



(3.33) 



(3.34) 
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Figure 5: Plot of polytropic k-essence 0;+ versus t with W = 0.2, B = —0.1, a 
-1.5 and (3 > 0. 

The pressure is given by 



a 



P 



1 — a V 1 — a 



-2 
l + o, 



{WX)'i^±J{WX)'i^ + ^ 

1 — a 



2 

1 + a 



B 



2a 
1 + a 



1 — a 

The corresponding EoS parameter is 

2 



(WX)— ± d{WX)^ + 



AB 
1 — a 



2a 
' 1 + a 



(3.35) 



a 



a — 1 



AB 



1 — a 



(WX) 



1 + a 

4,a ± 



(WX)"^ + 

1 — a 



-2 



(3.36) 



Using Eq. fl2.15p and after some manipulation, the EoS parameter takes the form 

5(2PF/32t2(/3-l) 



a 



1 + a 
I 2a 



[(a- l)(Vr/32t2(/3-i))^ _2^5]2' 



(3.37) 



Figure 5 shows the same behavior of EoS parameter as discussed for n = 2. 
However, this remains in the region, —1 < u+ < —1/3, at different points. In this 
case, the EoS parameter becomes polytropic k-essence EoS parameter as a < 0. It 
is concluded that this kind of parameter shows the accelerating expansion of the 
universe. 



3.0.5 Example 5: n = —1 

For this example, A and a satisfy the relation 

2a 



A 



A 



or a 



2(1 + A)' 



2a + 1 

Consequently, Eq. (13.51) becomes 

(1 + v4)pi+" - (iyX)-^p-(i+") -5 = 0, 



(3.38) 



(3.39) 
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which yields the solution 



P 



B±\l + A{1 + A){WXy 
2(1 + A) 



2a 



or 



P 



The pressure is 



B{2a + l) ^ 2a + l 



2a + I 



{wxy 



l + a 
2a 



(3.40) 



(3.41) 



2a 



P 



2a + 1 



B(2„ + 1)^2«+2,/b. + ^(„,;,)- 

2a + 1 



l + a 
2a 



B 



2 2 V 2a + 1 

The corresponding EoS parameter yields 



l + a 

2a 



2a 



2a + 1 



B 



2a + 1 



l + a 
2a 



-1 



(3.42) 



(3.43) 



Inserting the value of kinetic term X, the positive and negative roots of EoS param- 
eter u become 



2a 



UJ4 



52(W^/32t2(/3-l))^ 



2a + 1 52(2a + l)(iy/32t2(/3-i))^+2^' 



2a 



2a + 1 



+ B' 



l + a 
2a 



(3.44) 
(3.45) 



The left graph in Figure 6 shows the cosmological evolution of EoS parameter 
a;+ with respect to time. It represents a quintessence region of the universe for 
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higher values of /3 while for smaller values, it becomes constant and behaves like 
cosmo logical constant. In the right graph, polytropic k-essence EoS parameter is 
plotted. It represents initially a matter dominated universe for t < 1.2 for smaller 
values of /3. Then cj_ crosses the phantom divide line from quintessence to phantom 
region and becomes constant in the phantom region (cj_ < — 1). 

4 Other Forms of k-essence Cosmology 

There are two other possibilities of solvable k-essence cosmology [12] such as: 

• General k-essence cosmology: K{X, (p) = V{(f))F{X) 

• Quintessence k-essence cosmology: K{X, 0) = V{(f)) + F{X) 

These forms yield the same equations (13.51) and (13.61) for p and for parameter n. It 
is mentioned here that the above forms show the same behavior as in the case of 
purely kinetic k-essence cosmology. 

Finally, we would like to mention the other class of solvable k-essence models 
namely integrable k-essence models [20]. As an example, let us consider the following 
so-called E//£;-model 

Kxx = 2K^ + XK + ai^), (4.1) 

where a = a{(j)). It is nothing but the Painleve II (P-II) equation which is known 
as integrable. Let the solution of this equation we define as i^' = K{X] a) and here 
we assume that a = constant. Then the equation (4.1) has the following particular 
solutions (see [20] and references therein) 

K = K{X;1.5) = iIj -{2iIj'^ + X)-\ (4.2) 

K ^ KiX;l) = -^, (4.3) 

1 3 X'^ 

K = KiX- 2) = — - , , (4.4) 
T-^/x^ ON 6X2(X3+10) 

K = K(X-^-\- ^ I 6-^'(-^' + 10) 9X^(X^ + 40) 

^ ' ' X X6 + 20X3 - 80 X9 + 60X6 + 11200' 
K = K{X;0.5e) = -eij (4.7) 

and so on. Here 

^ = (Inx)x, X(X,0) = C,i(t>)Ati-2"'/'X) + C2i(l>)Bti-2-'/'X), (4.8) 

where Ci = Ci{(j)) and Ai{x),Bi{x) are Airy functions. Each of these solutions 
determine different type of cosmology. For example, the solution (4.4) corresponds 
to the EoS: 

2(2 -X3) ^ ^ 

6(X<> - lOX' - 8) ,, 
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Hence we get the EoS parameter 



1 



(4.11) 



3 X6-10X3-8' 



Inserting the value of kinetic term (2.15), this equation takes the form 



1 



3206/36t6(/3-l) 



(4.12) 




,12^12^12(/3-l) _ 8O06/36t6{/3-l) _ 512- 



If /9 = 1 = 00; then u ~ —0.27. This shows that the solution (4.4) represents a 
universe which is nearly to meet the quintessence phase of the accelerating universe. 
Similarly, we can calculate the EoS and its parameter for the other solutions from 
the system (4.2)-(4.7) and as well as for the other integrable k-essence models from 
|20] . Finally we note that it is interesting to extend results of this work to the 
periodical and quasi-periodical gas/fluid models (see e.g. Ref. [21]) as well as to the 
knot universe models [22j . 

5 Conclusion 

It has always been interesting to explore models to study the accelerating expansion 
of the universe. This paper is devoted to the study of some solvable k-essence 
cosmologies using MCG. In particular, we have discussed the purely kinetic k-essence 
model by taking the energy density and pressure of MCG model. We have found 
a nonlinear differential equation for the energy density which depends upon the 
parameter n. We have constructed EoS parameter for different values of n by using 
power law form of potential. Also, the values of constants involved in the model 
such as. A, B, a and /3 are constrained in such a way that they give the compatible 
behavior of our expanding universe with acceleration. The graphical results can be 
summarized as follows: 

• Example 1 is taken as n = which further leads to two cases. The case 
a = gives accelerating expansion of the universe dominated by cosmological 
constant for /3 > 0, while /3 < indicates a universe when DE is dominated 
over the matter to accelerate the universe. The second case, A = —1, indicaress 
the dynamics of the universe underlying in quintessence phase. 

• The second example leads to the generalized Chaplygin gas model for n = 1. 
For this value of parameter, initially our universe shows quintessence behav- 
ior of DE dominated era. After that, it converges to u = —1 representing 
cosmological constant. 

• For n = 2 and (3 < 0, the EoS parameter shows quintessence behavior of the 



• The case n = 0.5, using k-essence polytropic EoS parameter, gives the same 
dynamics of the universe as discussed for the parameter n = 2. 



universe. 
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• For n = —1, UJ+ leads to a cosmological constant dominated universe for 
smaller values of (3 whereas the k-essence polytropic EoS parameter U- repre- 
sents a quintom behavior and becomes constant in the phantom region. 

It is concluded that this relationship provides consistent results about acceler- 
ating expansion of the universe [13]. We would like to mention here that the same 
behavior is observed for the general k-essence and quintessence cosmologies as dis- 
cussed for puerly kinetic k-essence cosmology by using MCG model. Finally, we have 
given some integrable k-essence cosmologies K{X,a). These models show different 
phases of the accelerating universe depending upon its parameters. 
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